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ABSTRACT 

We consider the linear axisymmetric stability of a differentially rotating colli- 
sionless plasma in the presence of a weak magnetic field; we restrict our analysis 
to wavelengths much larger than the proton Larmor radius. This is the kinetic 
version of the magnetorotational instability explored extensively as mechanism 
for magnetic field amplification and angular momentum transport in accretion 
disks. The kinetic calculation is appropriate for hot accretion flows onto com- 
pact objects and for the growth of very weak magnetic fields, where the collisional 
mean free path is larger than the wavelength of the unstable modes. We show 
that the kinetic instability criterion is the same as in MHD, namely that the an- 
gular velocity decrease outwards. However, nearly every mode has a linear kinetic 
growth rate that differs from its MHD counterpart. The kinetic growth rates also 
depend explicitly on /3, i.e., on the ratio of the gas pressure to the pressure of the 
seed magnetic field. For (3 ~ 1 the kinetic growth rates are similar to the MHD 
growth rates while for (3 3> 1 they differ significantly. For (3 ^> 1, the fastest 
growing mode has a growth rate ~ for a Keplerian disk, larger than its 

MHD counterpart; there are also many modes whose growth rates are negligible, 
^ (3~ l l 2 VL <C Vt. We provide a detailed physical interpretation of these results 
and show that gas pressure forces, rather than just magnetic forces, are central 
to the behavior of the magnetorotational instability in a collisionless plasma. We 
also discuss the astrophysical implications of our analysis. 
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1. Introduction 

Balbus & Hawley (1991; BH91) showed that differentially rotating accretion disks are 
linearly unstable in the presence of a weak magnetic field (see Balbus & Hawley 1998 for 
a review; BH98). This instability, known as the "magnetorotational instability" (MRI), is 
local and extremely powerful, with a growth rate comparable to the rotation frequency of the 
disk. MHD turbulence resulting from the MRI is the most promising source of the efficient 
angular momentum transport needed in astrophysical accretion flows (e.g., Hawley, Gammie, 
& Balbus 1995; Armitage 1998; Hawley 2000; Stone & Pringle 2001). The MRI may also be 
important for the dynamo generation of galactic and stellar magnetic fields. 

In this paper we present a linear analysis of the MRI in a collisionless plasma using 
kinetic theory. The kinetic calculation is appropriate whenever the wavelength of the un- 
stable modes is shorter than the collisional mean free path. This regime is astrophysically 
interesting for several reasons: 

(1) In MHD, the most unstable mode of the MRI has a wavelength A as va/Q, where 
va = B/y/Airp is the Alfven speed and f2 is the rotation frequency of the disk. Thus, for a 
very weak magnetic field, the fastest growing mode has a very short wavelength, less than 
the collisional mean free path in many cases. A kinetic treatment is therefore required to 
determine whether the MRI can amplify very weak fields (e.g., the "first" magnetic fields 
generated at high redshift by a Biermann battery or analogous mechanism). 

(2) Radiatively inefficient accretion flows onto compact objects provide a useful framework 
for interpreting observations of low-luminosity X-ray binaries and active galactic nuclei (see, 
e.g., Ichimaru 1977; Rees et al. 1982 [the ion torus model]; Narayan & Yi 1995 [ADAFs]; for 
a review see Narayan et al. 1998 or Quataert 2001). In such models, the accreting gas is a 
hot two-temperature plasma in which the proton temperature (~ 10 12 K near a black hole) 
is much larger than the electron temperature (~ 10 9 — 10 11 K). In order to maintain such 
a two-temperature configuration, the accretion flow must be effectively collisionless in the 
sense that the timescale for electrons and protons to exchange energy by Coulomb collisions 
is longer than the inflow time of gas in the accretion disk. In principle, a kinetic treatment 
of the accretion flow structure, rather than a fluid treatment, is therefore necessary. The 
calculations described in this paper represent a first step towards understanding the physics 
of angular momentum transport and the structure of the accretion flow using kinetic theory. 

Our analysis is restricted to wavelengths much larger than the proton Larmor radius 
and frequencies below the proton cyclotron frequency. To motivate why kinetic effects can 
be important even on these "large" scales, consider a uniform medium threaded by a weak 
magnetic field (/3 > 1, where f3 is the ratio of the gas pressure to the magnetic pressure). 
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There are three long-wavelength waves in such a system: (1) the sound wave, (2) the Alfven 
wave, and (3) the slow magnetosonic wave. It is well known that the sound wave and the 
slow wave are very different in a collisionless plasma than in collisional plasmas described 
by MHD (e.g., Barnes 1966). We shall see that the same is true for the MRI. 

This paper is organized as follows. In the next section we describe our basic equations 
and assumptions (§2). In §3 we discuss linear waves in a collisionless plasma, emphasizing 
an important difference between MHD and kinetic theory that is useful for understanding 
the kinetic MRI results. In §4 we numerically solve the kinetic MRI dispersion relation and 
discuss its physical interpretation. We also show that a generalization of Balbus & Hawley's 
(1992; BH92) "spring" model of the MRI captures the main results of the kinetic calculation. 
Finally, in §5 we summarize our results and discuss their astrophysical implications. 



2. Basic Equations 

In the limit that all fluctuations of interest have wavelengths much larger than the 
proton Larmor radius and frequencies much less than the proton cyclotron frequency, a 
collisionless plasma can be described by the following fluid equations (e.g., Kulsrud 1983; 
Snyder, Hammett, & Dorland 1997): 

^ + V-( P V) = 0, (1) 

<9V (V x B) x B 
P^ + P(V-V)V = ^ -J. V-P + F g , (2) 

^Vx(VxB), (3) 

and 

P = p±I + (p\\ - p±)bb, (4) 

where p is the mass density, V is the fluid velocity, B is the magnetic field vector, F g is the 
force due to gravity, b = B/|B| is a unit vector in the direction of the magnetic field, and I 
is the unit tensor. Equations (jl[)-0 are identical to the basic equations of (collisional) MHD 
except that the pressure, P, is a tensor that is generally different perpendicular (p±) and 
parallel (j>\\) to the background magnetic field (e.g., the temperature need not be isotropic 
in a collisionless plasma). Formally, the pressure in equation (Q) should contain a sum over 
all particle species in the plasma (electrons, protons, and ions). In what follows, however, 
we consider a single fluid model in which only one contribution to the pressure response 
is included. In practice, the ions dominate the dynamics under consideration and so the 



-4- 



pressure can be interpreted as the ion pressure. This is particular true for hot accretion 
flows in which T p 3> T e . 

In a collisionless plasma, the parallel and perpendicular pressures satisfy separate "equa- 
tions of state" given by (e.g., Chew, Goldberger, & Low 1956): 

" B sfe) = - v ' (6 ^-^ v - 6 (5) 

and 




B 2 dt 



-V-(bgn)-2g ± V-b, (6) 



where d/dt = d/dt + V • V is the Lagrangian derivative, and q± and q» represent the flow 
of heat due to the motion of particles along magnetic field lines. Note that although there 
is no heat flow perpendicular to the magnetic field due to the very small proton Larmor 
radius, the perpendicular pressure/temperature can change due to heat transport along the 
magnetic field and so q± ^ 0. If one neglects the heat flux terms, equations (H) and (|) 
reduce to "double adiabatic theory" (Chew et al. 1956). Equation @ then describes the 
invariance of the average magnetic moment of the plasma, /i oc T±/ B oc p±/(pB), where 
T± is the perpendicular temperature. And equation (|6]) describes adiabatic parallel pressure 
changes due to the expansion or contraction of fluid elements (Kulsrud 1983). 

Equations ([[])- (§|) can be rigorously derived by expanding the Vlasov equation in the long 
wavelength, low frequency limit, and taking velocity moments (e.g., Kulsrud 1983). They 
face, however, the usual problem that the heat fluxes q± and qu depend on the third moments 
of the particle distribution function and so additional equations are needed to "close" the 
moment hierarchy. In Kulsrud's kinetic MHD one avoids this closure problem by solving 
the drift-kinetic equation, which is the low-frequency, long-wavelength, limit of the Vlasov 
equation (see Kulsrud 1983). By taking moments of the resulting distribution function one 
calculates p± and p\\ for use in equation (Q). For linear problems this approach is not too 
difficult and is the one employed here (see eqs. ]TJ| and [jTSJ below). For nonlinear problems, 



however, it is much more involved. Snyder et al. (1997) developed fluid approximations 
for q± and q\\ that model kinetic effects such as Landau damping and phase mixing. In this 
approach one solves equations @ and @ instead of solving for the full distribution function. 
For nonlinear problems, this is computationally more efficient and is a possible way of using 
MHD codes to extend the linear results of this paper to the nonlinear regime.^ 



5 Although we use the full drift-kinetic equation to calculate p± and pn, we have also found that the 
closures in Snyder et al. (1997) provide an excellent approximation for the linear problems considered here. 
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2.1. Linear Perturbations 

We assume that the background (unperturbed) plasma is described by a non-relativistic 
Maxwellian distribution function with equal parallel and perpendicular pressures (tempera- 
tures). Although the equilibrium pressure is assumed to be isotropic, the perturbed pressure 
will not be, which is a crucial difference between the kinetic and MHD problems. We take 
the plasma to be differentially rotating, but otherwise uniform (e.g., we neglect temperature 
and density gradients). Thus the velocity is given by V = Vo + 5v where Vo = RQ<j) an d 
Q(R) is the rotation rate. We consider a weak (subthermal) magnetic field with vertical 
[B z = B sin9) and azimuthal (B^ = B cos8) components, where 9 = ta.n~ l [B z / B^] is the 
angle between the magnetic field vector and the direction and B Q is the magnitude of the 
seed field. In a differentially rotating plasma, a finite Br leads to a time-dependent B^, 
which greatly complicates the kinetic analysis (unlike in MHD, where a time-dependent B^ 
can be accounted for; e.g., BH91); we therefore set B R = 0. Finally, we consider fluctuations 
of the form exp[— iut + ik ■ x], with k = k R R + k z z, i.e., axisymmetric modes; we also restrict 
our analysis to local perturbations for which |k|i? 1. Writing p = p + Sp, B = B + 5B, 
P± — Po + and p\\ = Po + 8p\\, and working in cylindrical coordinates, the linearized 
versions of equations (0)-© become: 

u5p = p k ■ <5v (7) 

- iujpoSvR - p 2Q5v (t> = (B Z SB Z + B^SB^) + lkzBz5Br _ ik R Sp ± (8) 

47T 47T 



K 



2 ik 7 BJB, 



- iujpodv^ + p 5v R — = — ik z sin 9 cos 9 \Sp\\ - 5p±\ (9) 

— iup 5v z = — z ^ — - — ik z sin 2 98p\\ + cos 2 95p± (10) 

u5B r = —k z B z 5v R (11) 
ik B dVL 

ljSB^ = -k z B z 8v<h — — — -5v R + B^k ■ 5v (12) 

uj dmR 

u5B z = k R B z 5v R , (13) 

where k 2 = 4f2 2 + dQ 2 /d In R is the epicyclic frequency. Equations (|7|)-(p^D are very similar 
to the analogous equations in BH91 except that we do not impose incompressibility and 
the pressure response is anisotropic. In particular, note that, even though we consider 
axisymmetric modes, there is a pressure force in the 0- momentum equation (eq. [[J) because 
the perturbed pressure is anisotropic (i.e., • (V • P) = ik z P z ^ = ik z (5p\\ — Sp±) sin^cos^). 

To complete our system of equations we need expressions for Sp± and 8p». These can 
be obtained by taking the second moments of the linearized and Fourier transformed drift- 
kinetic equation and are given by (e.g., eqs. 23-25 of Snyder et al. 1997) 
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Po Po \ B 

and 



^ (14) 



S P , SP, D (S_P_6B\ 



Po Po \Po Bq 

where |B| = B + 5B, 5B = b • <5B is the parallel magnetic field perturbation, and 



D 1 = l-R(£), D 2 



(16) 



Note that the second terms on the right hand side of equations ( |i~4D and (|T5|) are the per- 
pendicular and parallel temperature perturbations. In equation ([TBI) , R(£) = 1 + £Z(£) is 
the plasma response function, 



is the plasma dispersion function (e.g., Stix 1992), and £ = u> / (\/2c \k\\ |), where k\\ = b ■ k 



is the wavevector along the magnetic field and Co = yT/m is the isothermal sound speed of 
the particles (we have absorbed Boltzmann's constant into T so as to not cause confusion 
with the wavevector). 

Because equations ([L4]) and (|l|) are rather different from the MHD equation of state it 
is worthwhile discussing their physical interpretation. Consider first fluctuations for which 
£ > 1, in which case Z(£) « -£~ x - 0.5£~ 3 - 0.75£- 5 , R(£) w -0.5£~ 2 - 0.75£~ 4 , D 1 « 1, 
and -D 2 ~ 2. Equation ([TJ|) thus reduces to 

^^ + ^ (18) 

PO PO B 



and equation (p~5|) reduces to 

^3^-2^. (19) 
Po 5 

These are the linearized double adiabatic equations (eqs. || and || with q± = q\\ = 0). 
Not surprisingly, the adiabatic limit requires uj 3> fcyco, i.e., that the fluctuation timescale is 
much less than the time it takes particles to stream across the wavelength of the mode. 

In the opposite limit, £ <§C 1, Z{£) « iy^r ~~ 2£ and i?(£) « 1 + i^C, — 2£ 2 , so that 
Di ~ D 2 « —i^fn^. Equations fllj]) and (|l|) thus reduce to 

^JjL-i^if) (20) 

Po PO \Bq) 
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and 



(21) 



Po Po \Po B J 

These correspond to nearly isothermal fluctuations: the temperature perturbation is smaller 
than its "natural" value by a factor 1. This is because uj <C k\\Co, i.e., particles stream 

across a wavelength on a timescale much less than the fluctuation timescale and efficiently 
wipe out temperature gradients. Equations ( p0|) and ([H]) are the appropriate limit for the 
MRI. This is because the MRI has ^ k z v& and (3 1, and so £ ^ /J^ 1 / 2 ^ 1. 



2.2. The Dispersion Relation 



To obtain the dispersion relation, we eliminate all non-velocity variables from the mo- 
mentum equations. We first calculate 5B = cos 95 + sin 95 B z using 5B^, from equation 
(|T2|) and 5B Z from equation (|13|): 



SB 2n k ' Sv -a n k z 5v <t> 
— — = cos 9 sins' cos # - 

Bq UJ UJ 



d£l . a Q k z 5v n . ., 
— - sin cos h sin 9 



:'"'/>' :_2nk R 5v R 



dlnR a; 2 ' ' ^ 

Substituting equations (^) and fl22|) into equations ( |i~4"l) and (pT5|) then yields and 5p\\ as 
functions of 5v, which we substitute into the perturbed momentum equations (eqs. |p|-|T0|]). 
Using the perturbed induction equations (eqs. [p| -[p^|) we can also eliminate £B from the 
momentum equations in favor of 8v. This yields the following versions of the momentum 
equations: 

= 8v R [u 2 - k 2 v 2 Az - k 2 R (v 2 M + cl) + ik R k z VAzVM + ik R k z c 2 — sin 9 cos 9 f ^ 

(J (lliifi uj a In _rt 

k 2 R c^Di\ + 8v^ k R k z VA z VA<j> — 2iflu + k R k z c 2 ] D 1 sin 6* cos 6* 

k R k z (cl + v 2 M ) + A; R fc 2 Co cos 2 0£>i , (23) 







8v. 







r in 2 uj 
2?T 



dlnR 



+ v A 0V A zk R k z + k z k R c 2 D 1 sin 9 cos # 



,iM (jD2 + jDl ) sin 2^ cos 2^ d(J 



5va 



uj dlnR 

+ <5t> 2 k 2 v A<t> v Az — k 2 c^D 2 sin 3 9 cos 6* + k 2 c^Di sin cos 3 # 



.k 2 c 2 



5v R [-k R k z {v 2 M + cl)+lk 2VMVAz ^ < ' ~- 2 " 



(24) 

dn 



,, , -fc z fcflc;^icos 2 + i-^Dicos d 0siii0- r 
a; am it uj am it 



;^ J D 2 sin 3 ecos0^,l+^ 



a; 



dlnR 1 



k z VAcj>VAz — k 2 z (?^D 2 sin 3 9 cos + k^c^Di cos 3 6* sin 6 1 



2„2; 



+ 5v z uj 2 - k 2 z {v 2 A s + eg) - A; 2 c^ 2 sin 4 - fc;cg.Di cos 4 



(25) 
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where 2 

v\ z = -^- and v 2 M = ^- (26) 

are the Alfven speeds associated with the vertical and azimuthal fields, respectively. Equa- 
tions (|23|)-(Po]) define a matrix equation of the form A5v = 0. Setting det(A) = gives the 
dispersion relation. We have not found it particularly illuminating to write out the entire 
dispersion relation, nor have we made much progress solving it analytically, so instead we 
proceed to discuss its numerical solution. We will also present a simple model problem that 
captures the essential physics of the kinetic MRI. 

The MHD dispersion relation for the MRI, including the effects of compressibility, can 
be obtained from equations fl2"B|)-([2l)D by setting D\ = D 2 = 0. Equations (|T^) and ( [To] ) 
show that, for D x = D 2 = 0, 5p± = 5p\\ = Spc^, i.e., the perturbations are isothermal and 



the perturbed pressure is isotropic. Our basic linear perturbation equations (eqs. @-[|3 
reduce to their MHD analogues in this limit. In particular, note that in MHD the MRI is 
independent of whether the perturbations are adiabatic or isothermal; this is because it is an 
incompressible instability so the precise form of the sound speed is irrelevant for /3 ^> 1 (e.g., 
BH91). Thus the key simplification to the kinetic equations obtained by setting Di = D 2 = 
is that the perturbed pressure becomes isotropic, as it is in MHD. 



3. Linear Waves in Double Adiabatic Theory 

Before considering the full kinetic MRI problem, it is instructive to consider the sim- 
pler problem of linear waves in a uniform medium. In particular, we show that the slow 
magnetosonic wave is very different in kinetic theory than in MHD. This is important for 
understanding the kinetic MRI because the slow wave, along with the Alfven wave, is central 
to the dynamics of the MRI. We use double adiabatic theory throughout this section. Al- 
though double adiabatic theory does not include collisionless damping, which is quite strong 
for the slow mode and would alter some of the results in this section, it does show the signif- 
icant differences introduced by the anisotropic pressure in a collisionless plasma. Since our 
interpretation of the kinetic MRI in §4 focuses on the importance of this anisotropic pressure, 
it is useful to see its implications first in a simpler problem. In §4 and the Appendix we show 
that the qualitative conclusions drawn in this section carry over to the full kinetic analysis. 

Double adiabatic theory in a uniform medium corresponds to setting Q = k = and 
£ ^> 1 in equations (|2^ ) - (|25|) , in which case D\ = 1 and D 2 = 2. Without loss of generality 
we can take = so that cos 6 = and sin# = 1. To make contact with standard notation, 
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we also write k r = k± } k z = k\\, and va z — Va- The dispersion relation is then given by 



z ; 2 2 
UJ — K\\V 



lj 2 - k v A - 2ki4){u A - 3k 2 c 2 ) - kikicl = D A D MS = 0. (27) 



The analogous MHD dispersion relation is 

u 2 - k\v 2 A ] [(u 2 - k 2 v\ - klv 2 )(uj 2 - k 2 v 2 ) - klk^t] = 0, (28) 

where v 2 s = 7Cq is the adiabatic sound speed and 7 = 5/3 is the adiabatic index. 

Equation shows that, as in MHD, the double adiabatic dispersion relation factors 
into two parts: an Alfven wave branch (Da = 0) and a magnetosonic branch (Dms — 0). The 
Alfven wave in double adiabatic theory is identical to that in MHD, while the magnetosonic 
waves are different - this is because the "adiabatic index" in a collisionless plasma is different 
for motions perpendicular and parallel to the magnetic field. Motion along the field is 
one-dimensional and corresponds to 7 = 3 (hence the 3k 2 term in eq. |p7[] ) while motion 
perpendicular to the field is two-dimensional and corresponds to 7 = 2 (hence the 1k\ term 
in eq. P7[| ). By contrast, in MHD, the pressure is isotropic and 7 = 5/3. 

Figure 1 shows a plot of the dispersion relation of the fast and slow magnetosonic 
waves in MHD (dotted lines) and in double adiabatic theory (solid lines), taking (3 = 100. 
The fast wave, which is essentially a sound wave, is qualitatively similar in the two cases 
(the quantitative differences are due to the different 7's). The slow wave, however, is quite 
different. In MHD, the dispersion relation of the slow wave is degenerate with that of the 
Alfven wave for /? > 1, namely u = k\\VA- Except for k± = 0, this is not true in double 
adiabatic theory. The frequency of the slow wave depends on the sound speed; in fact, for 
fci / 0, the primary restoring force for the slow wave in double adiabatic theory is gas 
pressure, not magnetic forces. 

This result can be understood as follows. In MHD, the properties of the > 1 slow 
wave can be calculated by explicitly imposing incompressibility, V • 5v oc 5p ~ 0. This 
additional constraint (incompressibility) replaces the equation of state to determine the 
pressure (the Boussinesq approximation). In a collisionless plasma, this cannot happen 
because the pressure response is different parallel and perpendicular to the magnetic field, 
i.e., there are two equations of state (one for p± and one for py). Both equations of state 
cannot be replaced by the single requirement that the fluctuations be incompressible. More 
physically, a k± ~ kn slow wave in MHD has 5p ~ 5B 2 /8tt ~ B 5B±/4tt, i.e., the gas 
pressure, magnetic pressure, and magnetic tension forces are all comparable. Equivalently, 
Sp/po ~ (3~ 1 5B/Bq ^ 5B/B (for /3 ^> 1). In double adiabatic theory, however, a parallel 
magnetic field perturbation 5B/B Q induces a pressure perturbation Sp± t \\/po of comparable 
magnitude (see eqs. [JTSj and fl9[| ). This means that the pressure forces are much larger than 
the magnetic forces [5p ~ (3 5B 2 /8ir 3> SB 2 /8ir) and dominate the dynamics of the wave. 
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The exception to these arguments is if the pressure perturbation vanishes, i.e., 5p± = 
8p» = 0. Alfveh waves and the k± — > limit of the slow magnetosonic wave are the only 
waves in MHD that have 5p = (they also have 8B = 0). As a result, these pressure-free 
waves are the only incompressible fluctuations in double adiabatic theory. For all other 
waves, and in particular for slow waves with k± ^ 0, pressure is the dominant restoring force 
in a (3 ^> 1 plasma and so the frequencies depend strongly on the sound speed (Fig. 1). 

The results in this section are relevant to the MRI because the MRI is an incompressible 
instability with \uo\ <C kcQ. Although pressure forces generally lead to a small modification 
of the MRI in MHD, they will be substantially more important in the kinetic analysis (just 
as for the slow magnetosonic wave in this section). 



4. The Kinetic MRI 

As noted in §2.3, the general kinetic MRI dispersion relation appears to be analytically 
intractable. In this section we present its numerical solution and physical interpretation. 
As a check on our numerical calculations we have confirmed that our results reproduce the 
kinetic dispersion relation for the Alfven wave and the slow and fast magnetosonic waves 
when Q = (including the collisionless damping rates) .0 We also reproduce the MRI in MHD 
when the kinetic terms are dropped (this requires setting D 1 = D 2 = in eqs. p3|-p5|). 



Figures 2-4 show the results of numerically solving the kinetic MRI dispersion relation, 
assuming a Keplerian disk for which f2 oc R~ 3 ^ 2 and k — Q. The figures show the kinetic 
growth rate of the MRI for different values of f3 z = 87cp /B 2 , for different magnetic field 
geometries (defined by B<p/B z ), and for different wavevectors (ku and k z ). The corresponding 
MHD results are shown for comparison by the dotted lines. It is important to note that in 
MHD the MRI growth rate is essentially independent of j3 and B^/Bz] by contrast, Figures 
2-4 show that the kinetic results depend sensitively on both of these parameters. 

Figures 2-4 show that, although the growth rates can be very different, the region 
of instability in wavevector space is the same in kinetic theory and MHD. To understand 
this result, it is sufficient to consider the u — > limit of the kinetic equations, since this 
determines the transition between stable and unstable modes. Setting u = implies that 
£ = u / (y/2k\\Co) = as well. From equations fll4D -(|l6|), it then follows that 5p±/po = 
$P\\/Po = dp/Po- Physically, as £ — >• 0, there is more and more time for particles moving 



6 We compared our results to the linear kinetic code described in Quataert (1998) and to the analytic 
results in Barnes (1966) and Foote & Kulsrud (1979). 
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along magnetic field lines to efficiently transport heat. This leads to nearly isothermal 
fluctuations in which the pressure perturbation is isotropic and is set only by the density 
perturbation. As discussed in §2.3, the kinetic equations reduce to the MHD equations in 
this limit. This is an important result because it shows that the MHD instability criterion, 
namely dQ 2 /dR < 0, applies to the kinetic problem as well. Moreover, the set of unstable 
modes is the same in MHD and kinetic theory, as is seen explicitly in Figures 2-4. 

Perhaps the three most striking results of the kinetic calculation shown in Figures 2-4 
are: (1) The kinetic growth rates depend sensitively on (3. For (3 ^> 1 they differ significantly 
from the MHD growth rates while for (3 ~ 1 they are similar (see Fig. 3). (2) For = 0, or 
for sufficiently large the kinetic growth rates are smaller than their MHD counterparts, 
particularly at large (3 (e.g., Fig. 3a). (3) For ^ 0, the kinetic growth rates can be larger 
than their MHD counterparts (e.g., Fig. 2 and Fig. 4b). Moreover, for (3 ^> 1, the fastest 
growing mode is at k z v Az *C where there is negligible growth in MHD (Fig. 4b). 

To understand the kinetic MRI results, we have found it useful to consider the equa- 
tions that describe the displacement of a fluid element from its equilibrium circular orbit. 
BH92 and BH98 showed that, for the special case of a vertical magnetic field and vertical 
wavevector, the radial and azimuthal components of the MHD momentum equation can be 
written in terms of the radial and azimuthal fluid displacements, £r and as 

~W - 2 ^ = - [d^R + {kzVAz) ) (29) 

^ + 2^ = -(k z v Az )%. (30) 
As discussed by BH92 and BH98, equations (|29"D and ( p0|) are identical to the equations 



describing two orbiting point masses connected by a spring of spring constant k 2 v 2 Az (in 
MHD, magnetic tension plays the role of the spring). This suggests the following physical 
interpretation of the MRI in MHD (BH92). For a rotation profile with dVt 2 jdR < (unstable 
to the MRI), a fluid element at radius R — 5R is rotating slightly faster than a fluid element 
at radius R. The "spring" pulls backwards on this inner fluid element, removing its angular 
momentum and forcing it to move to a yet smaller radius. Similarly, a fluid element at radius 
R + 5R is rotating slightly slower than a fluid element at radius R and so the "spring" pulls 
forward on this fluid element, giving it angular momentum and forcing it to move to a yet 
larger radius. This simple physical picture captures the essence of the MRI in MHD. 

A useful toy model that provides additional insight into the physics of the MRI, both 
in MHD and kinetic theory, is given by the following equations for the fluid displacement 

~w ~ 2n ^r = ~ [ihR + Kr (31) 
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Equations ( pTj) and (|32|) describe the displacement of rotating fluid elements coupled by an 
anisotropic "spring," for which the spring constant is different in the azimuthal (K^) and 
radial (Kr) directions (this is clearly no longer a real spring!). The unstable root in the 
dispersion relation associated with equations (|31|) and fl32|) is given by 



2 k + Kr + K<j> 1 
U = 2 2 



n\2 / dfi 

K 4> + K R + k J - 4i^ AT fl 1 



1/2 

(33) 



For /Tr = = (k z VAz) 2 , equation ([33|) gives the k,R = dispersion relation of the MRI in 
MHD (this is shown by the dotted line in Fig. 4b). It is also straightforward to show that, 
for Kr > K^, the growth rates in equation (^) are smaller than the MHD growth rates (i.e., 
those with Kr = K$) and for > Kr the growth rates in equation ( |33"D are larger than 
the MHD growth rates. For example, for = 0, equation (|33p gives u = for any Kr and 
so there is no instability. On the other hand, the ^> Q ^> Kr solution of equation ( |33"1) 



is \u\ ~ y | dVL 2 /d In R\. For a Keplerian disk this gives |u;| = v^^? which is larger than the 
growth rate of the fastest growing mode in MHD = 30/4). 

These results can be understood physically by noting that it is ultimately the presence 
of an azimuthal restoring force, rather than a radial restoring force, that is destabilizing in 
the MRI. This is because it is the azimuthal force that removes angular momentum from 
an inwardly displaced fluid element and adds it to an outwardly displaced fluid element. 
By contrast, the radial force is stabilizing because it attempts to "pin" the fluid element 
to its equilibrium position. Thus K^ > Kr leads to faster growth because it enhances the 
destabilizing azimuthal force relative to the stabilizing radial force (and vice- versa for Kr > 
K^). For the remainder of this section we explain how thermal pressure in a collisionless 
plasma plays the role of the anisotropic "spring" in the above toy model. This will account 
for the behavior of the kinetic MRI seen in Figures 2-4. 

Because they are restricted to Ur = and — 0, equations (|29| ) and ( pOf ) do not 
include the effects of gas pressure or magnetic pressure (both of which vanish in this special 
case). To understand the kinetic MRI we need to include these restoring forces using the 
radial and azimuthal momentum equations (eqs. || and 0, respectively). This yields the 
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following equations for the fluid displacement^ 

-W- 2Q ^ = -{dhui + {k ^ ) J^H^^J' (34) 

^ + 2^ = -(k z v Az )% - ik z {^j^j ^6cose. (35) 

In equations (|34[) and (|35|) we have simply rewritten the pressure gradients from equations 
(|) and @; in the Appendix we calculate these explicitly in terms of the fluid displacement. 
It is worth noting again that there is a pressure force in the 0-momentum equation (eq. |35[] ) 



because the perturbed pressure is anisotropic. In MHD, 6p\\ = 5p± and so this term vanishes. 

Following the arguments in §3 and Figure 1 we expect that the pressure gradients in 
equations (B3J) and fl3"5] ) will be much more important in kinetic theory than in MHD. In the 
Appendix we calculate the magnitude of these pressure forces and confirm this hypothesis. 
We use these results below to present a physical interpretation of the kinetic MRI, focusing 
on two important special cases: (1) = 0, ^ 0, for which the kinetic growth rates are 
smaller than their MHD counterparts (e.g., Figs. 3a & 4a), and (2) fcjj = 0, B^ ^ 0, for 
which the kinetic growth rates are larger than the MHD growth rates (e.g., Fig. 4b). 

Consider first the special case of = and Ur ^ (e.g., Fig. 3a & 4a). In this case a 
displaced fluid element feels a restoring force in the radial direction due to gas and magnetic 
pressure; there is, however, no analogous pressure gradient in the direction (only magnetic 
tension). This corresponds to Kr > in the toy model of equation (33); the growth rates 



should therefore be suppressed with respect to the k R = growth rates. The presence of a 
stabilizing radial pressure gradient provides a physical explanation for why the MHD growth 
rates decrease with increasing kR (see, e.g., the dotted line in Fig. 3a). Moreover, in the 
Appendix we show that the pressure gradient in kinetic theory is larger than in MHD by a 
factor of ~ /3 1 / 2 . The kinetic growth rates should therefore be even smaller than the MHD 
growth rates, with stronger suppression at larger (3. This is precisely what is seen in the 
kinetic calculation; e.g., Figures 3a and 4a show the = growth rate for different (3. 

Consider now the special case of k B = 0, but B^ ^ (e.g., Fig. 4b). In this case the 
radial pressure force vanishes, but there is an azimuthal pressure force due to the anisotropic 



pressure. As suggested by the toy model in equation (|33|) this azimuthal pressure force 



which is not present in MHD, is destabilizing because it removes angular momentum from 



7 Strictly speaking, equations (|3^) and ( |35| ) should have an additional term on the right hand side given 
by k z VAzVA0(k ■ 0- For (3 3> 1, this term is negligible because the MRI is nearly incompressible and so we 
do not consider it further. 
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an inwardly displaced fluid element and adds it to an outwardly displaced fluid element 
(just as the azimuthal component of magnetic tension does). Moreover, for ~ B z the 
destabilizing pressure force is larger than the destabilizing magnetic tension force by a factor 
of ~ (3 l l 2 (see the Appendix). This explains why the Ur = 0, B^ ^ growth rates are larger 
than their MHD counterparts, and why the growth rates increase with increasing f3 (see, e.g., 
Fig. 4b). It also explains why the growth can be rapid even at k z VA z ^ when magnetic 
tension (which drives the MRI in MHD) is very weak. In fact, for f3 ^> 1 and k z VA z "C fi, 
the forces in the kinetic MRI are arranged as follows: azimuthal pressure ^> Coriolis ^> 
magnetic tension. We therefore expect the growth rates to approach the ^> Q ^> K R 
limit of equation (^), namely |u;| ~ y/3Cl. As shown in Figure 4b, the fastest growing modes 
do approach this maximal growth rate. 

Although the above interpretation focuses on two special cases, the results in Figures 
2-4 can be readily understood as a competition between the stabilizing radial pressure force 
and the destabilizing azimuthal pressure force. The importance of gas pressure, rather than 
magnetic forces, also explains why the kinetic results depend sensitively on (3. 

5. Summary and Discussion 

In this paper we have presented a linear axisymmetric calculation of the magnetorota- 
tional instability (MRI) in a collisionless plasma. Our analysis is restricted to wavelengths 
much larger than the proton Larmor radius, frequencies below the proton cyclotron frequency, 
and "seed" magnetic fields with no radial component (Br = 0). The MRI is believed to give 
rise to MHD turbulence and efficient angular momentum transport in astrophysical accretion 
flows, and may also be important for the dynamo generation of galactic and stellar magnetic 
fields (e.g., BH98). Our kinetic calculation, rather than an MHD calculation, is appropriate 
whenever the collisional mean free path of the protons exceeds the wavelength of the MRI. 

The instability criterion for the kinetic MRI is the same as in MHD, namely that the 
angular velocity decrease outwards. The set of unstable modes is also the same in kinetic 
theory and MHD. However, nearly every mode has a linear kinetic growth rate that differs 
from its MHD counterpart. For example, the fastest growing mode in kinetic theory has 
a growth rate ~ \^3£l for a Keplerian disk, which is larger than its MHD counterpart by 
a factor of 4\/3/3 ~ 2.3.0 More generally, the kinetic growth rates can be either larger or 
smaller than the MHD growth rates, depending on the orientation of the magnetic field and 
the wavevector of the mode (Fig. 2). The kinetic growth rates also depend explicitly on /3, 



This rapid growth is obtained only for Urvaz <§C ^, k z VAz *C fi, /? 1, and ^ B z (see Figs. 3 & 4). 
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i.e., on the ratio of the gas pressure to the magnetic pressure. For (3 3> 1 the kinetic results 
differ significantly from the MHD results while for (3 ~ 1 they are similar (see Fig. 3). 

We have argued that the kinetic MRI can be understood by considering the force due 
to pressure gradients in a high (3 collisionless plasma. In MHD, pressure leads to a relatively 
minor modification of the MRI. In kinetic theory, however, the pressure forces are ~ /3 1 / 2 
times larger than in MHD and are therefore dynamically much more important (see §3 and 
the Appendix). Moreover, in kinetic theory there is an azimuthal pressure force even for 
axisymmetric modes (so long as ^ 0; see eqs. |L0| and p5| ). This is because the pressure 



response is anisotropic in a collisionless plasma: it is different along and perpendicular 
to the local magnetic field. This azimuthal pressure force, which is not present in MHD, is 
destabilizing because it removes angular momentum from an inwardly displaced fluid element 
and adds it to an outwardly displaced fluid element (just as the azimuthal component of 
magnetic tension does in MHD). The destabilizing pressure force explains why the kinetic 
growth rates of the MRI can be larger than the MHD growth rates (e.g., Fig. 4b). 

The importance of gas pressure shows that the character of the MRI is somewhat dif- 
ferent in a collisionless plasma than in a collisional plasma described by MHD. The crucial 
function of the magnetic field is to enforce an anisotropic pressure response, rather than to 
directly destabilize the plasma via magnetic tension. The importance of pressure gradients 
also explains why the kinetic results depend sensitively on (3. For (3 ~ 1 pressure forces are 
comparable to magnetic forces, and the kinetic growth rates are not that different from the 
MHD growth rates, while for (3 ^> 1 pressure forces dominate over magnetic forces and the 
kinetic results differ substantially from the MHD results (e.g., Fig. 3). 

BH92 showed that the MRI in MHD could be understood using a simple model in 
which magnetic tension acts like a spring coupling different fluid elements in the plasma. 
We have presented a generalization of BH92's "spring" model that captures many of the 
results of the full kinetic MRI calculation (see eqs. @]-^3[). In this model the radial and 
azimuthal "spring constants" are different; physically, this corresponds to the anisotropic 
pressure response in a collisionless plasma. 

To conclude, we briefly discuss the astrophysical implications of our results, focusing 
on the two applications mentioned in the introduction: (1) the amplification of weak fields 
generated by a Biermann battery or analogous mechanism, and (2) hot two-temperature 
accretion flows onto compact objects. 

(1) For a very weak magnetic field MHD predicts that the fastest growing mode of the MRI 
has a very small wavelength « va/Q oc B. This will be less than the collisional mean 
free path in many cases. Our kinetic analysis shows that there is rapid growth of the MRI 
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even in this limit. This is encouraging for the hypothesis that the MRI contributes to the 
dynamo amplification of very weak magnetic fields, e.g. the generation of galactic fields from 
a cosmological seed field. To further assess this question, however, it is necessary to extend 
our analysis to include finite Larmor radius effects. In particular, the "battery" generation 
of magnetic fields is limited by self-induction to field strengths such that the proton Larmor 
radius is comparable to the size of the system (e.g., Balbus 1993). Finite Larmor radius effects 
will always be important on these scales, particularly since the wavelengths of unstable MRI 
modes are then much less than the proton Larmor radius. 

(2) In radiatively inefficient accretion flows onto compact objects, which have been applied 
extensively to low-luminosity accreting sources (e.g., Narayan et al. 1998), the inflowing gas 
is a hot two-temperature plasma in which the proton temperature is much larger than the 
electron temperature. In order to maintain T p 3> T e , the timescale for electrons and protons 
to exchange energy by Coulomb collisions must be longer than the inflow time of the gas. 
This requires a sufficiently small accretion rate, M ^ o?Medd (e.g., Rees et al. 1982), 
where Medd is the Eddington accretion rate and a is the dimensionless Shakura-Sunyaev 
viscosity parameter. Since the timescale for proton-electron collisions to modify the proton 
distribution function is comparable to the proton-electron energy exchange timescale, the 
proton dynamics is effectively collisionless for any two-temperature radiatively inefficient 
accretion flow;[] the kinetic calculation presented in this paper is therefore appropriate for 
describing angular momentum transport by the MRI in such mo dels .p°] 

It is, however, difficult to apply our linear calculations to the nonlinear saturated state 
expected in the accretion flow. Nonetheless it is worth noting that there are rapidly growing 
modes in a collisionless plasma even for (3 ^> 1 so weak fields can be efficiently amplified. 
Moreover, MHD simulations find saturation at (3 ~ 1 — 100 with a predominantly toroidal 
field (e.g., BH98; Stone & Pringle 2001). For this magnetic field configuration, the linear 
kinetic growth rates of the MRI are not that different from their MHD counterparts (if 
anything, they may be somewhat larger; e.g., Figs. 3b & 4b). While this suggests that 
the saturated turbulence may be qualitatively similar in kinetic theory and MHD, there will 
undoubtedly be quantitative differences. In addition, the fact that the fastest growing modes 
occur at somewhat different wavenumbers could change the nonlinear results. Perhaps more 
importantly, collisionless damping of the sound wave and the slow magnetosonic wave is very 
strong and operates on all scales in a collisionless plasma, while strong damping in MHD is 



Proton-electron collisions are more important than proton-proton collisions because T p S> T e . 

10 By contrast, geometrically thin accretion disks (e.g., Shakura & Sunyaev 1973) are much cooler and 
denser; MHD accurately describes the dynamics of thin disks so long as the gas is sufficiently ionized. 
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restricted to very small scales. This may alter the nonlinear behavior of the MRI. Numerical 
simulations that address these issues would be extremely interesting. 

Our results may also have implications for understanding particle heating in radiatively 
inefficient accretion flows. The radiative efficiency of such models is set by the amount of 
electron heating in the plasma. This depends on how the energy in MHD turbulence is 
dissipated (e.g., via a turbulent cascade, reconnection, etc.). The prominent role of pressure 
fluctuations in the kinetic MRI suggests that the resulting turbulence may couple better to 
slow waves (which have a pressure perturbation) than Alfven waves (which do not). Slow 
waves primarily heat the protons in the collisionless plasmas of interest (e.g., Quataert 1998; 
Blackman 1999) while an Alfvenic cascade may lead to significant electron heating if (3 ^ 10 
(e.g., Gruzinov 1998; Quataert & Gruzinov 1999). Kinetic simulations of the MRI should be 
able to assess the relative importance of slow wave and Alfven wave excitation. 

We thank Steve Balbus, Steve Cowley, Barrett Rogers, Alex Schekochihin, and Anatoly 
Spitkovsky for useful discussions. GH was supported in part by the U.S. Department of 
Energy under Contract DE-AC02-76CH03073. EQ was supported in part by NASA grant 
NAG5-12043. 
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A. Calculation of the Pressure Forces 

In this Appendix we calculate the radial and azimuthal pressure forces in equations ( p4|) 
and ( |55| ) in terms of the fluid displacements £r and These are used in our interpretation of 
the kinetic MRI results in §4. We restrict our analysis to the two important limits highlighted 
in §4: (1) = 0; k R and (2) y^0;k R = 0. 



A.l. B^ = 0;k R ^0 
In this case there is a radial pressure force given by (eq. 



'A 
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SB 2 
8np 



5p± 
Po 



-ik 
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S P 2 . SB ( 2 n 2 

— c + — { v Az + D lCo 
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(Al) 



where we have used Sp±/p = Sp/p + Di5B/B from equation ( |14|) in the second equality. 
We now rewrite all of the terms in F R in terms of £r, the radial displacement. For B^ = 0, 
SB = 5B Z . The radial component of the induction equation (eq. [fLl]] ) thus yields 



SB _ k R 5v R 
B u 



-ik R C 



R- 



(A2) 



To calculate Sp/po = {k z 5v z + k R 5v R )/uj in terms of ^ R alone we need to find 5v z as a 
function of <5t> R . To do this note that the z-component of the momentum equation (eq. ) 
implies 

k ¥p\\ 



k z Sv z 



Lup 



(A3) 



Since Sp/po — 5B/B = k z 5v z /u, equation ([15]) gives 5p\\ as a function of both Sv z and Sv R . 
Substituting this into equation ( |A3| ) we solve for 5v z in terms of 5v R and thus find 



r 2 k 2 



5p ., / 

PO V w - C k z( l + ^2) 



(A4) 



Substituting equations ( |A2|) and (|A4|) into equation (|A1|) , and assuming /3 ^> 1 so that 
|cj 2 | k 2 cl, yields 



+n z 1- 



a; 



fcM.(l + 2£> 2 ) 2 . 



(A5) 



The MHD limit of equation ( |A"5| ) can be obtained by setting Di = D 2 = (see §2.3). In 
this case F R = —£ R k R v Az [l + \u\ 2 /(k 2 v 2 Az )} ~ — £Rk R v Az . Consider instead double adiabatic 
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theory, for which Di = 1 and F* 2 = 2. In this case Fr ~ — £r^r: c o; this is larger than the MHD 
pressure force by a factor of ~ (3. Finally, for the full kinetic problem we need to evaluate 
Di and D 2 using equation fllTf). Since the MRI has \w\ & k z VA we can take £ <C 1 so long as 
/3 ^> 1. In this case D\ « _D 2 ~ —iV^£ ~ —iu}/k z Co, so that Fr ~ —k 2 R ^ R c^{—iuj/k z CQ). To 
estimate the magnitude of Fr, note that |cj| ~ A; z t>^ in MHD, in which case Fr ~ — ^Rk 2 R v a 2 c . 
This is ~ /3 1//2 times larger than the pressure force in MHD. This large radial pressure gradient 
suppresses the growth rates of the MRI, as seen in Figures 2-4. 



A.2. B^^ 0; k R = 



In this case there is an azimuthal pressure force given by (eq. ||35| 



-ik 7 sin 9 cos 9 
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where we have used equations ([14]) and (|T5|) to eliminate <5p^ and foy. For > 1 the MRI 
is nearly incompressible and Sp/p <C 5B/B .f^ We therefore neglect the Sp/p term in 
equation (|A"6|). Using = d^/dt — ^RdVt/d\n.R one can rewrite equation (|22"D for 55 in 
terms of Again neglecting 5p/po relative to the other terms, this yields 



— = — ik z sin 9 cos 
Bo 



Substituting equation (|A7|) into equation ( |A6|) yields 



F = -^A; 2 c 2 sin 2 9 cos 2 9(D 2 + D x 



(AT) 



(A8) 



In MHD, F<£ = 0, and magnetic tension, which ~ — k 2 z v 2 Az ^ (see eq. |35fl), plays the 
destabilizing role. By contrast, in kinetic theory the azimuthal pressure force is given by 
F<$> ~ — ^ sin 2 9 cos 2 9k 2 cl(—iuj/k z c ). For ~ £> z , so that sm9 ~ cos^ ~ 1, this is larger 
than the destabilizing azimuthal tension force by a factor of ~ (3 1 / 2 . This large destabilizing 
azimuthal pressure force enhances the growth rates of the MRI, as seen in Figures 2-4. 



n The calculation in Appendix A.l shows this explicitly: eqs. [A4] and ]A2{ imply that &p/ pa ~ D28B/ Bq 
where D 2 ~ (3~ l/2 < 1. The same qualitative result holds for the different geometry considered here. 
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Fig. 1. — The dispersion relation for the magnetosonic waves in MHD (dotted lines) and in 
double adiabatic theory (solid lines), taking (3 = 100. In MHD, the slow wave dispersion 
relation is identical to that of the Alfven wave (uj = k\\v a) while this is only true for k± <C k\\ 
in double adiabatic theory (see the text for an explanation). 
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Fig. 2. — The kinetic growth rates of the MRI for f3 z = 8irpo/Bl = 10 4 and for different 
geometries of the seed magnetic field. The MHD results, which are independent of B^, are 
shown for comparison (dotted line). The vertical wavenumber is taken to be h z v\ z j^l = 
y!5/16, which is the fast growing mode in MHD. 
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Fig. 3. — The kinetic growth rates of the MRI for varying f3 z (solid lines). The MHD results, 
which are nearly independent of f3 z , are shown for comparison (dotted line). The vertical 
wavenumber is taken to be k z VA z /Q = ^/ 15/16, which is the fast growing mode in MHD. 
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Fig. 4. — The kinetic growth rates of the MRI as a function of k z for different (3 Z (solid lines). 
The corresponding MHD results are shown by the dotted line. 



